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We present exact solutions, mainly analytical, for the two-site double-exchange-Holstein model, that 
allow us to draw a complete picture of the role of both phonon and spin quantum fluctuations in 
determining the short-range correlations in the manganites. We provide analytical solutions of the 
model for arbitrary electron-phonon coupling and phonon frequency, for S = 1/2 and for the classical 
spin limit S — oo, and compare these results with numerical diagonalization of the realistic S — 3/2 
case. The comparison reveals that the realistic case S = 3/2 is not well described by the classical spin 
limit, which is often used in literature. On the other hand, the phonon fluctuations, parametrized 
by the phonon frequency wo, stabilize ferromagnetic phases with respect to the adiabatic limit. We 
also provide a complete analysis on the polaron crossover in this model. 
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I. INTRODUCTION 

It is known from the 1950s that the double exchange 
mechanism |l|-|3| is at the basis of the magnetic proper- 
ties of the manganese perovskites Ri_ x A x Mn03 (where 
R is a rare earth element (e.g., La), and A is a divalent 
element like Sr or Ca). In this compounds, the <i-levels of 
each Mn ion host 4 — x electrons. Three of them occupy 
the three low-lying t2 g levels, with aligned spins due to 
the Hund's rule. These electrons are basically localized, 
and only the 1 — x electrons in the e g level contribute to 
the transport properties. The discovery of the so-called 
"colossal" magnetoresistance has originated an enor- 
mous revival of studies on these compounds, both on the 
theoretical and the experimental side. The systematic 
experimental investigations of the last few years have 
underlined some weakness in the previous understand- 
ing, unveiling a surprisingly rich phase diagram where a 
lot of competing phases are stabilized by varying doping, 
temperature and chemical nature of the dopants. One of 
the most relevant new theoretical trends is the suggestion 
that the transport properties cannot be fully understood 
on the basis of the double exchange alone, and that the 
interplay of this mechanism and a significant electron- 
phonon (e-ph) interaction leading to a Jahn- Teller effect 
is the key for the explanation of these properties || . The 
preminent role of e-ph effects has also been firmly estab- 
lished experimentally by various groups and techniques 

The complex entanglement between charge, orbital 
and lattice degrees of freedom represents a hard theo- 
retical challenge, that is far from being solved. Many 
approximate solutions and numerical results have been 



proposed, but the complexity of the phase diagram has 
naturally forced various authors to many uncontrolled 
simplifications. In this work we make a step back, and 
focus on a system for which analytical exact results can be 
obtained. Namely, we solve the double exchange model 
for a single electron on two sites in the presence of a lo- 
cal Holstein e-ph coupling J7|. Since we are interested in 
the relevant physics determining the interplay between 
lattice and spin quantum fluctuations, we consider the 
simple Holstein model, instead of a more involved Jahn- 
Teller coupling. This choice does not imply a loss of 
generality since we are not discussing the role of orbital 
degrees of freedom. 

The two site system has been extensively studied as 
the minimal system able to capture the key features of 
polaron formation from the point of view of ground state 
[pHilf as well as spectral properties |l^,|l3[. Recently 
also the interplay between e-ph and e-e correlations have 
been studied semi-analytically within the same model 
fll4f . However, as far as magnetic properties are con- 
cerned, it is quite obvious that a two-site system does not 
allow for long-range order and phase transitions. Never- 
theless, it shows short-range (nearest neighbors) corre- 
lations, that give substantial indications on the actual 
long-range properties of the system, at least in strong 
coupling. In the context of the models for the mangan- 
ites it is in fact believed that the finite-size effects play 
a little role pi . In the following we will define " first or- 
der transitions" , the level crossing between phases with 
different symmetry, and "second order transitions" the 
continuous transitions. We will discuss the relevance of 
our results to large systems in the following. 

Due to the simplicity of the model, we can give com- 



plete exact phase diagrams without approximations. One 
of our main results is a complete characterization of the 
role of the quantum fluctuations of the core t 2g spins. In 
a microscopic model of the manganites, the core spins 
have S = 3/2, and this value is usually thought to be 
large enough to get rid of quantum fluctuations and treat 
them as classical variables. We will explicitly test this as- 
sumption by comparing the limiting cases S — oo (classi- 
cal spins) and the extreme quantum case S — 1/2, where 
the effect of quantum fluctuations is maximum, with the 
realistic S = 3/2 case. 

Analogously, we will discuss the role of lattice quan- 
tum fluctuations, releasing the adiabatic approximation 
on the phonon degrees of freedom. The role of quantum 
phonon fluctuations is not trivial, as already known for 
e-ph models alone 16 18). For S = oo and S = 1/2 we 
give analytical exact solutions of the model, exploiting an 
exact analytical solution of the two-site Holstein model. 
For S — 3/2 we use standard exact diagonalization to 
solve the model. Also in this case no approximation is 
introduced and all the regimes are accessible. 

A similar study has been reported in Ref. 119], where 
the two-site double exchange model for classical spins is 
solved by perturbation theory around a variational ref- 
erence state obtained by a Lang Firsov canonical trans- 
formation. Our work overcomes some limitations of Ref. 
[19|, namely the classical spin limit. Contrary to Ref. 
[19 1, we are also able to explore the adiabatic regime 
luq <C t, well beyond the region in which the Lang Firsov 
result is a good reference state. 

The paper is organized as follows: In Section O we 
introduce the model and the methods for our analytical 



solutions; In Section III we present the phase diagram of 



the model and discuss the role of quantum fluctuations; 
due to the complexity of the phase diagram the discussion 
is divided in three subsections: in the first the effect of 
the magnetic degrees of freedom on the polaron crossover 
is considered; in the second we discuss the effect of the 
e-ph interaction on the magnetic phase diagram of the 
model, and in the third subsection the full phase diagram 
is presented. 



In Section IV we discuss the relevance of our results 



for larger size systems and for the experimental scenario. 
Finally we give concluding remarks in Section [V]. 



II. METHODS OF SOLUTION 

We consider the Holstcin-Double Exchange model on 
a two-site cluster for a single electron: 

H = -t ^(c T 1(T C2 :(J + c\ a d, a ) + 
—Jh 2_^ Ui ' ^ + ^i^ 1 ' ^ 2 + 

8=1,2 

— g(m — n 2 )(a + a') + ujoa^a, (1) 



where Si (i = 1,2) is a local spin associated to the lo- 
calized t 2g electrons on each site, Ci,<r(c\ a ) destroys (cre- 
ates) an electron of spin a on site i, n» = X)<x C I o- c i>" 
is the number operator on each site, <Xj = c ia a a pCip is 
the spin operator on each site (a are the Pauli matri- 
ces), a (a 1 ) is the destruction (creation) operator for a 
lattice distortion that couples to the difference of den- 
sity between the two sites. The lattice displacement X is 
given by X — y/h/2muio(a + a'). We could have started 
from a standard Holstein model with a phonon mode per 
each site, coupled to the local density. It is in fact easy 
to show that, in this case, the symmetric combination 
of the two phonon modes A = 1/V2(ai + a 2 ) couples 
to the total density, giving rise to a trivial term, and 
the only term left is the one we introduced, where the 
phonon mode may be written in terms of the local ones 
as a = l/-y/2(ai — a 2 )- 

We explicitly consider, besides the hopping between 
the two sites and the Hund's rule term (Jh) that couples 
ferromagnetically the conduction electrons to the local- 
ized ones, an antiferromagnetic superexchange term (Ji) 
between the core electrons. Jh is always taken to be the 
largest energy scale, consistently with the physics of the 
manganites. This latter term, even though J\ is signif- 
icantly smaller than Jh, has crucial importance on the 
magnetic properties of the manganites pO] . We also con- 
sider a Holstein coupling (g) between the electron and a 
dispersionless mode of frequency uj . 

In this paper, we present the exact solution for the 
model (|l|). In particular for the classical spin case S = oo 
and the extreme quantum case S — 1/2 we provide ana- 
lytical solutions for arbitrary values of both the electron- 
phonon coupling and of the phonon frequency ujq, ex- 
ploiting an exact solution of the two-site Holstein model 
reported in Appendix [Aj. In the S = 3/2 case, a numeri- 
cally exact solution by means of exact diagonalization is 
instead presented. In the following subsections, we de- 
scribe the analytical solutions for S = 1/2 and S = oo. 

A. Exact Solution for S — oo 

Following Ref. j||, in the classical spins case we can 
write 



H(8) = Ji cos 6 - H H (6) + E 



Hund 



(2) 



where Enund is the contribution of the Hund's term to 
the total energy and Hh is the hamiltonian of a 2-site 
Holstein model in which the hopping t is replaced by 
t = tcos(6/2). The canting angle 8, that measures the 
relative orientation of the core spins fully characterizes 
the magnetic arrangement. If 9 — the spins are aligned 
and a ferromagnetic (FE) state is found, whereas for 
= it, an antiferromagnetic state (AF) is found. Inter- 
mediate values of 9 describe canted (CA) states. There- 



fore, the solution of the classical spins two-sites Holstein 
double-exchange model can be obtained by minimizing 
on 8, once the eigenvalues of Hh are known, as shown in 
Appendix 0. For Jh S> 3\ the extremal condition then 
gives 



sin(0/2)(-2Jicos(0/2) 



tdE 



H 



2 dt 



)=0. 



(3) 



which shows that the ferromagnetic state 9 = is always 
an extrema of E(9) pi] . Then the transition from FE to 
CA state is of "second order". The critical coupling J\ 
for this transition is given by the vanishing of the term 
in parentheses in Eq. (0) 



•n 



Tpkin 



(4) 



where E H ln = tdEn/dt is the kinetic energy of the Hol- 
stein model, i.e., the kinetic energy of the system with 
J\ = (Notice that E^ is a negative quantity). The 
effect of el-ph interaction on the FE — > CA transition is 
therefore the substitution t — v —Ekm- 

Now let us consider the "first order" FE — v AF tran- 
sition. We have to compare the FE and AF energies 
obtained by Eq. (||) respectively with 6 = and 6 = n. 
The critical coupling is given by 



Jl = 



E H (0)-E H (t) 



(5) 



where Ejj(0) = —g 2 /uiQ is the energy of the atomic Hol- 
stein model. 



B. Exact Solution for S = 1/2 

In this case, neglecting for the moment the phonon de- 
grees of freedom, the electronic Hilbert space (including 
the core spins) is in principle made by 16 states, that 
reduce to 8 if the symmetry for inversion of all the spins 
is considered. As shown in Appendix M, this problem 
can be simplified, and the largest subspace to deal with 
is a 3x 3 sector, but the remaining problem is still not 
trivial if we switch on the coupling with the phonons. 
Fortunately, in the limit J# S> t, a further simplification 
occurs (also shown in Appendix |B|) , leading to the possi- 
bility to express the eigenvalues of the model in terms of 
the two-site Holstein model. The details of the solution 
are reported in Appendix |B|. In such a way, we can char- 
acterize the condition for the only possible transition, 
i.e., the transition from FE to AF ground state. The 
transition is discontinuous and can be obtained from the 
comparison of the energies of the different phases, that 
we compute in Appendix M. The critical coupling for the 
FE — > AF transition is then given by 



•/[ 



4{E H {t/2)-E H (t)) 



(6) 



where Eh is the energy of the 2-sites Holstein model. 
Contrary to the classical spin case quantum spin fluctu- 
ations allow for non zero effective hopping t = t/2 in the 
AF phase. 



III. RESULTS 

In this section, we present the phase diagrams of the 
model (ffl) , with a particular emphasis on the interplay be- 
tween the role of phonon and spin quantum fluctuations 
(measured respectively by ujo/t and the value of the "lo- 
cal" spin of the t^g electrons S) Due to the relatively large 
number of parameters that determine the phase diagram, 
we organize the discussion of the results in subsections: 
in the first subsection we discuss the polaron crossover 
in the different regimes, providing a unifying picture of 
the effect of both phonon and spin fluctuations, and of 
electron-spin correlation, that generalizes in a consistent 
way the conditions for the small polaron crossover in the 
simplest Holstein model [[16] Q. In the second subsec- 
tion, we discuss the magnetic transitions occurring in our 
model. Both the nature of the magnetic phases and the 
relation between the magnetic state and the occurrence 
of polaronic behavior are strongly dependent on the value 
of the spin S. Finally, in the last subsection, we present 
complete phase diagrams in the \-J\ plane, where the 
role of lattice and magnetic degrees of freedom is high- 
lighted. 



A. The polaron crossover 

Models with electron-phonon interaction quite gener- 
ally exhibit a polaronic ground state when the coupling 
strength is large enough. The transformation of the free 
electron into a small polaron is not a phase transition, but 
a continuous crossover. For the Holstein model and a sin- 
gle particle, it has been shown that the condition for the 
crossover significantly depends upon the ratio between 
the phonon frequency c^o, and the typical electronic en- 
ergy scale t. In the adiabatic regime ujq <C t, the crossover 
occurs for A = g 2 /u>ot ~ 1, whereas in the antiadiabatic 
regime luq 3> t, the crossover is controlled by the (purely 
phononic) variable a = g/(Vo |16H18|. As a matter of fact, 
the crossover coupling X po i is pushed to larger values of 
A as the frequency is increased. Moreover, the crossover 
becomes smoother and smoother as u)o/t gets larger. 

The conditions for a polaron crossover in the adiabatic 
and antiadiabatic regimes can be understood on basic 
physical grounds. In the adiabatic regime, the key con- 
dition is that a bound state can be formed. The condition 
A > 1 expresses this property, since it just implies that 
the polaron binding energy E po i = g 2 /ujq exceeds the ki- 
netic energy of a free electron ~ t. On the other hand, 
in the antiadiabatic limit, the electronic energy scale t is 



not the largest scale, and the polaron crossover is ruled 
by the condition a 2 > 1, that corresponds to the excita- 
tion of a significant number of phonons (or, equivalently, 
to a sizeable lattice distortion) . The crossover conditions 
we have just described are based on simple, model in- 
dependent, physical insights, and are therefore expected 
to basically hold, with some marginal changes, also for 
more complicated models like the double-exchange model 
we are considering. 

It must be noted that, since the formation of a polaron 
is not a phase transition with an associated broken sym- 
metry, there is some ambiguity in determining a physi- 
cally sensible clear-cut criterion for the polaron crossover. 
In most previous studies, including that of Ref. |lS|| , the 
crossover line has been drawn as the locus of the points in 
which some relevant expectation values, like the electron- 
lattice correlation function l-/NJ2i( n iXi} or the average 
number of phonons in the ground state change their be- 
havior. This kind of characterization has no problem 
in the adiabatic regime, where the crossover is rather 
sharp, but it is more questionable in the antiadiabatic 
regime. In this work we use a much more definite cri- 
terion, that is based on a qualitative difference between 
polaronic and non-polaronic states. Namely, we study 
the (quantum) probability distribution function for the 
displacement operator P(X) = (0\X)(X\0), where |0) is 
the ground state wave function and \X) denotes the state 
with displacement X. In the adiabatic limit loq = 0, 
the phonon degrees of freedom are described by classical 
variables, and no quantum fluctuations are present. The 
solution of the model involves a minimization of the elec- 
tronic groundstate as a function of X. As a result, the 
probability distribution is a single (or a few) ^-function, 
centred at the values that minimize the energy. More ex- 
plicitly, if the system is not polaronic, a single value of X 
minimizes the energy, while in the polaronic regime, two 
different minima are obtained. The polaron crossover is 
then associated with the coupling value in which a single 
(^-function leaves place to two symmetric peaks. As soon 
as the quantum fluctuations of the lattice are restored by 
introducing a finite phonon frequency ujq, the 6- functions 
broaden, but the qualitative features do not change. The 
polaronic regime is characterized by a bimodal distribu- 
tion, while the non-polaronic state present a unimodal 
distribution. In the particular case of a single particle 
on two sites, the polaronic regime presents two symmet- 
ric peaks at X = ±Xq, and the non-polaronic state in 
non-distorted, so that P(X) is peaked at X = 0. 

Fig. YU shows the evolution of P(X) varying A in the 
two-site Holstein model for Wo/t = 0.1 (representative of 
the adiabatic regime) and u>o/t = 4 (representative of the 
antiadiabatic regime). In both cases a smooth crossover 
occurs between a quasi-free electron state (unimodal dis- 
tribution) and a polaronic state (bimodal distribution). 
The figure also clearly shows that the crossover in indeed 
much sharper in the adiabatic case than in the antia- 



diabatic one. Furthermore, P(X = 0) in the polaronic 
region rapidly vanishes soon after the crossover in the 
adiabatic case, while it stays finite in the antiadiabatic, 
despite it is a local minimum in both cases. 
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FIG. 1. Phonon Displacement Distribution Function P(X) 

for the two-site Holstein model for u>o/t = 0.1 and uio/t — 4. 

The various lines in the two panels correspond to different 

values of A. For uoo/t — 0.1 A ranges from to 1.4, whereas 

for tu /t — 4 < A < 4. 

It is worth to emphasize that that the monomodal to 
bimodal crossover of P(X) has also been reported as sig- 
nature to a crossover toward a "polaronic" state also in 
studies of the Holstein model in the thermodynamic limit 
using Dynamical Mean Field Theory J23], p3[. 



B. Magnetic Correlations 

In this section we discuss the behavior of the magnetic 
correlations in the two-site double exchange model. It 
should be clear that a such a small system can not un- 
dergo phase transitions, and that we are only able to 
describe short range correlations. We parametrize the 
magnetic correlations between the two sites by means of 
the scalar product (Si • S2) between the core spins. 

Since our interest in the model is motivated by the 
manganites, we will always assume that Jh is the largest 
energy scale. In the S — 00 case we let Jh go to infinity, 
and in the finite spin case, we take Jh = 10t. 

In the absence of electron-phonon coupling, and as- 
suming that Jh is the largest energy scale, the direct an- 
tiferromagnetic exchange between the core electrons de- 
termines the magnetic properties of the system. For zero 
and small J\ , the spins are ferromagnetically aligned due 
to the Hund's coupling. Increasing Ji, antiferromagnetic 
correlations tend to appear. 

The nature of the spin correlations depends crucially 
on the value of the spin S, since the latter rules the pos- 
sible values of (Si ■ S2). More explicitly, in the classical 
spin case, (Si ■ S2) = S 2 cos{9), where the canting angle 



8 between the spins is a continuous variable. The fer- 
romagnet continuously evolves into a canted state as J\ 
is enhanced. The canting angle asymptotically tends to 
7T, corresponding to the antiferromagnetic state, as J\ is 
enhanced. Panel (a) in Fig. displays the dependence of 
(Si • S2) on J1S 2 for the classical spin case for Jh = 00. 
In the quantum case the total spin is given by 

(Si • S 2 ) = i(S(S + 1) - ft (Si + 1) - S 2 (S 2 + 1)) (7) 

and assumes only a few values. For Si = S 2 = 1/2,5 = 
and 1 are the two only possible values. For Si = S 2 = 
3/2, we can have four values (S — 0, 1, 2, 3). It must be 
noted anyway, that the total spin operator S 2 does not 
commute with the Hamiltonian (ffl), so that the energy 
eigenstates have no reason to be eigenstates of S 2 . 

An inspection to the results in the absence of electron- 
phonon coupling shows indeed that, for S = 1/2, the 
magnetic state abruptly varies, for J\ = Ji(FE — AF), 
from a fully polarized ferromagnet (FE) to an antiferro- 
magnetic (AF) state, which is not fully polarized. The 
transition is a level crossing between two states with dif- 
ferent symmetry. The exact value of (Si • S2) in this state 
depends on both Jh and J\ , even if the dependence on 
J\ (for J\ > Ji(FE — AF)) is really weak, as it appears 
in panel (c) of Fig. || 

The S = 3/2 case is strictly analogous to S = 1/2, 
and shows the fully polarized FE state, followed by three 
different combinations of the AF states. Also in this case, 
the precise values of (Si • S2) in the three states depend 
on Jh and J\ , and the dependence on J\ is really weak in 
each region. Moreover, the state with the largest negative 
correlation is really close to the full antiferromagnet. We 
label the two intermediate spin phases as canted 1 (CA1) 
and canted 2 (CA2), and the "most antiferromagnetic" 
as antiferromagnetic (AF ) tout court. The dependence of 
(Si • S2) on J\ is shown in Fig. g (b). We notice that the 
scale of Ji associated with the change in the magnetic 
structure is consistent with the experimental estimates 
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FIG. 2. (Si ■ Si) for S = oo,3/2, 1/2 (from top to bottom) 
as a function of Ji. 



C. The phase diagram 

In this section we discuss the interplay between the 
magnetic properties and the e-ph coupling and finally 
determine the phase diagram of our model. We tune the 
relevance of lattice and magnetic degrees of freedom, by 
varying the strength of the electron-phonon coupling A, 
and of the antiferromagnetic coupling between the core 
spins J\. Then we draw various phase diagrams in the 
A — Ji plane. Each of the diagrams is characterized by 
the values of the spin S and of the phonon frequency 
wo, that parametrize the relevance of quantum spin and 
lattice fluctuations respectively. We consider the S = 00, 
S = 3/2 and S = 1/2 cases, and UJo/t — 0.1 (adiabatic 
regime) and 4 (antiadiabatic regime) . Finally, we always 
assume that the Hund's rule coupling Jh is the largest 
energy scale. In the classical case, we take Jh = 00, and 
in the quantum cases, we use Jh = 10i. 

We denote "first order transitions" (level crossings be- 
tween states with different symmetries) with full lines and 
"second order transitions" (crossovers between phases 
with the same symmetry) with dashed lines. 



1. The effect of Ji on the polaron crossover 



In section III A, we have briefly described the condi- 



tions ruling the polaron crossover in the Holstcin model. 
In the adiabatic regime LUo/t <C 1, the condition for a 
polaron ground state is that the polaron binding energy 
E P oi = g 2 /uq is larger than the free electron kinetic en- 
ergy, measured by t. It is quite natural to generalize this 
condition to the double-exchange model, at least when 



the polaron crossover occurs between two phases that 
share a common magnetic state. In this case, we can 
replace the bare hopping t by the "magnetically renor- 
malized" kinetic energy at A = (I). Thus the crossover 
condition is determined by the condition 



A - 9 ~ 1 

A mq — = — J- 

U) t 



(8) 



In the S = oo case the magnetic hopping is given by 



t = tcos(-), 



(9) 



where 9 is the canting angle between the t2 g spins plpj. 
In the quantum cases S = 1/2 and S = 3/2, one can 
view the canting angle as a quantized quantity, that can 
assume only a few discrete values. We anticipate that 
these are not the quantized values of the semiclassical 
approximation. 

Regardless the value of S and u>o/t, for small values of 
Ji the ground state is always ferromagnetic due to the 
Hund's rule. A crossover occurs between a ferromagnetic 
itinerant electron and a ferromagnetic polaron. Within 
this region, the magnetic hopping is fixed to the free value 
t = t, and does not depend on J\. As a result, the model 
is completely equivalent to a two-site Holstein model, and 
the crossover is associated with a vertical line in the A— J\ 
diagram, as shown in all the phase diagrams (Figs. [|]4|, 
BtlOJl)- The crossover value of A depends only on the 
ratio uio/t, and moves from the A ~ 1 in the extreme adi- 
abatic limit, to A ~ 1.2 for uio/t = 0.1, to a substantially 
larger value (A ~ 3.46) for coo/t = 4, where the condition 
for the polaron crossover is close to a 2 ~ 1 (that implies 
A~4). 

Increasing Ji, phases with antifcrromagnetic correla- 
tion between the core spins appear. The nature of this 
phases depends on the value of S, as shown in section 
[II B. We start from the quantum cases, that present 
sharp level crossings at A = 0, where t sharply jumps 
following the magnetic correlations shown in Fig. @. If 
we neglect the really weak dependence on J\ of (Si • S2) 
within a given magnetic phase, the polaron crossover is 
controlled by the condition (|J) , where t is the value cor- 
responding to the actual magnetic phase. 

The exact results obtained as described in section 
U confirm this expectation, and the polaron crossovers 
among phases with the same magnetic correlation are in 
fact delimited by vertical dashed lines in all the diagrams 
for S = 1/2 and S = 3/2 (Figs. Hll|). The value of 
the crossover coupling obviously changes in the differ- 
ent magnetic phases. The FE state is the one with the 
largest kinetic energy due to the double exchange mech- 
anism, so that the critical A is the highest in this phase, 
and it decreases by decreasing the value of the magnetic 
correlation according to Eq. (||) (see, e.g., Fig. |[||). 
Notice that the effect of the value of (Si • S2) on the 



crossover coupling is much more evident in the adiabatic 
limit (Fig. ||), where the competition between the po- 
laron energy and the kinetic energy rules the crossover, 
than in the antiadiabatic regime (Fig. |I|), where the elec- 
tronic kinetic energy is not the most relevant quantity. 
In the extreme antiadiabatic limit uio/t — > 00 this de- 
pendence must completely disappear, since the kinetic 
energy plays no role in the crossover. 

The formation of polaron does not only occur as 
smooth crossover between states with the same magnetic 
correlation. Indeed, in the S — 1/2 case (see Fig. ||), 
if we continuously increase J\, we have that, between 
the vertical dashed lines corresponding to the polaronic 
crossovers within the FE and the AF phases, a first- 
order transition (level crossing) occurs from a ferromag- 
netic non-polaronic state and an antiferromagnetic po- 
laron. The interplay between the localizing effect of both 
the e-ph and the antifcrromagnetic magnetic interaction 
strongly favors the AF polaronic state with respect to 
the competing phases. 

Similar level-crossings occur for all the finite-spin 
cases, with more involved details depending upon the 
value of S and loq. 

For example, in the antiadiabatic regime u>o/t — 4, the 
FE—AF polaron transition occurs only in a narrow range 
of parameters (see Fig. g, compared to the adiabatic case 
||). This is simply due to the fact that, in this regime, 
the crossover values for A in the FE and the AF are very 
close. In the extreme antiadiabatic limit u>o/t — > 00 this 
region would indeed vanish. 

In the richer S = 3/2 case, the CA1 and CA2 states 
intrude between the FE and the AF at weak e-ph cou- 
pling. In the adiabatic regime (see Fig. H), no polaron 
crossover occurs within the canted phases, and both these 
phases undergo a first-order transition to the AF polaron 
Only the FE and AF phases display the usual polaron 
crossover. In the antiadiabatic case, besides the afore- 
mentioned reduction of the regions in which the polaron 
formation becomes first-order, canted polaronic states 
are stabilized by the phonon quantum fluctuations. (The 
"critical" frequency above which canted polaronic states 
appear is u>o/t — 1). 

In the S — 00 case, where canted phases with a contin- 
uous canting angle are stable, the polaron crossover is not 
represented by a vertical line, since the kinetic energy is 
a continuous function of J\ with t = tcosO{J\), and 9( J\) 
is the value of the canting angle at A = 0. Again, all re- 
sults are consistent with the condition (M). In this case, 
the phonon fluctuations play a somewhat qualitative role. 
In the extreme adiabatic limit the C A state undergoes a 
level crossing to the AF polaronic state. The e-ph inter- 
action and the antiferromagnetic coupling Ji cooperate 
to stabilize the AF polaron without forming a canted po- 
laron. As soon as we introduce a finite, but small LUo/t, 
a tiny slice of a canted polaronic phase appears to bridge 
between the canted state and the antiferromagnetic po- 



laron. In the antiadiabatic regime, the huge quantum 
fluctuations strongly favor a canted polaronic state, and 
the antiferromagnetic polaron appears only for A > 10 
and J\ > 2. 

The above results show that the effect of the mag- 
netic correlations on the small polaron crossover is influ- 
enced by the spin quantum fluctuations. In particular, 
the S = 3/2 case, which is relevant to the manganites is 
not qualitatively similar to the classical spin case, that 
is usually considered for simplicity. Many features of the 
S = 3/2 case are in fact direct consequences of the quan- 
tum nature of the spins, and are similar to the simplest 
quantum case S = 1/2. 



2. The effect of e-ph coupling on magnetic transitions 

In this section we analyze how the various magnetic 
transitions described in Section III B are influenced by 
the e-ph interaction. In the previous section we have 
found analytical results for the transition from FE to 
AF states for Ju — > oo. The expressions, given by Eqs. 
(pf) and (p), can be recast in the common form 



lim(l-7)i 



E H (t+0.-i)t)-E H (t) 



J[S 2 = 



E H {AF)-E H (FE) 

(&I&2) FE _ (Si-S2) AF 



(10) 



where E H (AF) (E H (FE)) is the energy of the Hol- 
stein model for the antiferromagnetic (ferromagnetic) 
spin alignment, and (Si • S2) in the different magnetic 
phases is computed at g = 0. In particular, the evalua- 
tion of Eh for a given phase simply amounts to find the 
groundstate of the Holstein model where the bare hop- 
ping t is replaced by t. This relation also holds for the 
case of S = 3/2, once the proper value for the kinetic 
energy in the antiferromagnetic phase t = £/4 is used. 

Eq. (|10| ) results from the competition between the 
magnetic energy balance, controlled by J\, and the "po- 
laronic" energy, i.e., the energy resulting from the e-ph 
coupling. 

Since the AF phase has always a smaller hopping with 
respect to the FE phase, the first qualitative effect of 
the e-ph interaction is to lower the energy of the anti- 
ferromagnetic phase with respect to the ferromagnetic 
one, therefore favoring antiferromagnetism. The region 
of stability of the FE phase is always shrunk by increas- 
ing A. More generally, the e-ph interaction favors phases 
with smaller values of (Si • S2), so that the boundaries 
of the different magnetic phases are always marked by 
downward curves in the Ji-A plane. 

In the limit of small t 1 the energy difference in the 
numerator of (|10|) reduces to the pure kinetic energy of 
the Holstein model close to the atomic limit. Defining 
7 = £/£, we can write 

lim t ^o 0M7*) ~ -M*)) = 



= (1-7)* 



dE H {t) 
dt 



(1-7)* 
-25#»(l- 7 ). 



(11) 



In the small t limit the condition ( |l0| ) reduces then to 

J! = n(S)(-E k H m ), (12) 

where k(S) is a constant that contains the factor 1 — 
7, and depends on the value of the spin S and on the 
transition under consideration, such that, in the absence 
of e-ph interaction the condition is simply Jf = n(S)t. 
This latter relation is analogous to Eq. (g), since the 
effect of the e-ph interaction results in the substitution 
t — ► —Eft™, but, contrary to that, it is valid only for 
t->0. 

This result gives valuable informations about the role 
of the retardation effects in the e-ph coupling. In general 
terms, the interaction mediated by the phonons is in fact 
retarded, and becomes instantaneous only if ujo/t —>■ 00. 
In such a limit, as we have shown, the kinetic energy rules 
the magnetic transitions. As soon as the approximation 
t — > is released, the retardation effects imply that Eq. 
© must be used. Notice that E H {AF) - E H (FE) = 
_jr;kin _|_ ^ where A is the difference between the e-ph 
interaction energies in th two phases, and turns out to 
be always positive. The overall effect of the retarded e- 
ph interaction is therefore to reduce the stability of FE 
phases with respect to AF and CA phases. 

In the S — 3/2 case all the magnetic transitions are of 
"first order", as described in section 1IIB, so that sim- 
ilar arguments can be applied, and Eq. (nfl) allows to 
compute the transition coupling, once the energy of the 
Holstein model and the magnetic energy of the appropri- 
ate phases are known. 

The situation is different only for the "second order" 
transition between the FE and CA phases in the classi- 
cal spin case. In this case the continuity of the transition 
implies that the energy difference between the phases is 
infinitesimal, so that Ejj n = —tdEu{t)/dt rules the tran- 
sition not only for small t, but for arbitrary values of t, 
as shown by Eq. (|l|). This preliminary analysis suggests 
a really important difference between the classical spin 
limit S = 00 and the quantum S = 3/2 case. 

Now we can give some description and interpretation 
of the exact phase diagrams in light of the above analysis. 

In the adiabatic limit uio/t = 0, for finite value of S, 
and for A < A po ;, the energy of the e-ph model alone 
does not depend upon A, so that the relative stability 
of the various magnetic phases is in turn expected to 
be A-independent. The magnetic transitions are there- 
fore associated with horizontal lines in the \-J\ plane. 
If we introduce phonon fluctuations, the kinetic energy 
depends upon A also before the crossover, and the bound- 
ary between the magnetic states acquires a finite slope, 
that becomes larger and larger by increasing uio/t. This 



behavior can be easily seen comparing Fig. || with Fig. 
[| (or |5| with ||). The case of the classical spin variables, 
where the magnetization is a continuous variable, can be 
understood in similar terms. The crossover between the 
FE and the CA phase in fact closely follows the curve 
obtained using Eq. (Q). 

As discussed previously, the phonon fluctuations al- 
ways favor FE phases with respect to AF and CA 
phases, as it can be seen comparing the phase diagrams 
for ujo/t = 0.1 with the corresponding with uio/t — 4 (at 
the same value of A) . In we increase the quantum fluctua- 
tions of the phonons, by increasing the phonon frequency 
uiq, the retardation effects are decreased, so that the lo- 
calization of the electrons is made more difficult. An 
enhanced mobility of the electron results in an enhanced 
stability of the FE phases due to the double exchange 
mechanism. 
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IV. RELEVANCE FOR LARGER SYSTEMS 

In this section we discuss the relationship between the 
two-site model and larger systems. There are not indeed 
so many solutions of the double exchange model in the 
presence of e-ph coupling, and most of them are forced to 
consider the classical limit for both the phonons and the 
core spins. In the case of numerical calculation the main 
problem in dealing with quantum phonons and spins is 
the enlargement of the Hilbert space (that becomes in- 
finite for the case of phonons). This problems are even 
harder if the orbital degrees of freedom are taken into 
account. In Ref. Jig] , a numerical analysis of the double- 
exchange model in presence of a Jahn- Teller coupling 
is performed for the three-dimensional structure of the 
actual compounds, and for a number of electrons cor- 
responding to the stoichiometric LaMnC>3. The three- 
dimensional structure and the realistic hopping matrix 
elements allow for different antiferromagnetic arrange- 
ments of the spins besides the ferromagnetic state. In 
particular, realistic values of the parameters give rise to 
the A-type antiferromagnetism found in the experiments 
p5|,p0| . Of course our two-site system is not enough to 
distinguish among the various kind of antiferromagnetic 
orderings, but the phase diagram in Fig. 2(b) of Ref. |lf| 
and our phase diagram for classical spins and for small 
phonon frequencies (Fig. |7|) are rather similar. The shape 
of the curves separating the FE phase to the "non-FE" 
phases look very similar, and the crossover line between 
a FE with no distortions and a FE polaronic phase is 
very close to the one between the FE and the FE(JT) 
phase. The agreement is not only qualitative, but also 
quantitative as it can be checked by direct comparison. 

The capability of the two-site cluster to capture the 
physics of both the polaron crossover and the magnetic 
transitions in a more realistic system is consistent with 
previous claims (see, e.g. fig]) of the irrelevance of finite- 
size effects for these systems. 



The x — 0.5 system has been considered in Refs. 
p5| , p6[ . Also in this case, the phase diagram is quite 
similar to Fig. [7|: A first-order transition separates a FE 
state to x\ox\-FE states, and within the FE phase a ver- 
tical line denotes a crossover associated to small polaron 
formation (or Jahn- Teller effect). 

The two-site cluster gives therefore results really close 
to larger systems when the latter are available. This 
suggests that our results are generically representative of 
larger sizes, also for parameters that are presently un- 
accessible to numerical simulations. We emphasize that 
this agreement is not accidental and can be easily ra- 
tionalized. As mentioned in the introduction, the small 
polaron crossover is well represented by a two-site clus- 
ter. This is essentially due to the extreme short-range 
character of the polaronic state. The capability of the 
two-site cluster to describe the crossover from a delo- 
calized state and a polaronic state is well documented. 
It must be noticed, however, that our small cluster can 
not determine whether the polaronic and the delocalized 
states are Fermi liquids or not. 

As far as the magnetic properties are concerned, the 
two-site cluster has no long-range correlation, but only 
short-range correlations. Nevertheless the features re- 
lated to local interactions such as the ones in Hamilto- 
nian (PP) are well represented. In this spirit we expect 
the results of our calculations to be somehow related to 
a Dynamical Mean Field Theory, where the local quan- 
tum fluctuations are exactly taken into account, while 
the spatial correlations are frozen |27[. Similarly to the 
Dynamical Mean Field Theory, our exact solutions faith- 
fully reproduce the physics of the model if the local and 
short-range effects are dominant, as it is the case for the 
manganites. More specifically the one electron solution of 
the two site Holstein model (see Appendix [A|) is surpris- 
ingly similar to the one electron solution of the Holstein 
model for an infinite lattice within Dynamical Mean Field 
Theory. Looking at a given site, the role played by the 
"effective quantum medium" of Dynamical Mean Field 
Theory is here simply played by the other site. 



V. CONCLUSIONS 

In this work the two-site double exchange model for an 
electron coupled with phonons is solved exactly for an ex- 
tremely wide range of parameters and physical regimes. 
For 5 = 1/2 and S = oo we give an analytical exact so- 
lution for arbitrary e-ph coupling and phonon frequency. 
For S = 3/2 the solution is obtained through standard 
numerical techniques. The availability of these solutions 
allows us to study the effect of both phonon and spin 
quantum fluctuations, and of their mutual interplay. 

This study, though limited to the extreme small size of 
the two-site cluster, is shown to be a good description of 



larger systems, since the relevant physics involves local 
and short-ranged quantities. 

One of our main results is a complete characteriza- 
tion of the effect of the double exchange and of an an- 
tiferromagnctic coupling between the ti g spins on the 
small polaron crossover. In this regard, we give an an- 
alytical estimate for the crossover coupling, given by 
A„, 9 = g 2 /ajo~t ~ 1, where t is the kinetic energy renor- 
malized by the magnetic effects. 

From a complementary point of view, we considered in 
detail the effect of the e-ph interaction on the magnetic 
properties of the system. In this case, we give an analyti- 
cal condition for the magnetic transitions in the presence 
of a finite A, given by Eq. (|10|). This relation can be sim- 
plified in the limit t —t 0, where there is no retardation 
effect. The comparison between the general case and the 
atomic limit allows us to quantitatively describe the role 
of retardation in stabilizing AF (or CA) phases. 

A comparison of the realistic S = 3/2 case with the 
classical spin case S = oo shows that this latter approx- 
imation does not reproduce some qualitative features of 
the phase diagram. A proper study of the manganites 
should therefore take into account the quantum nature 
of the core spins. 
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APPENDIX A: THE EXACT SOLUTION OF THE 
TWO-SITE HOLSTEIN MODEL 

In this appendix we sketch the solution trough a con- 
tinued fraction expansion of the two site Holstein model. 
A continued fraction solution has been already reported 
in literature [g8[ for a related model in the field of quan- 
tum optics. Here we derive the continued fraction ex- 
pansion for eigenvalues and eigenvectors of the two site 
Holstein model using a different method. 

The Holstein model for an electron on two sites can be 
written in a pseudo-spin representation in terms of the 
Pauli matrices 



Hh = LUolaJa — gc z (a + a T ) — t<r x 



(AI) 



where 1 is the unity matrix. The hamiltonian can be 
diagonalized in the electron subspace using a transfor- 
mation introduced in Ref. |p| 



\/2 



I (-) a a 
-1 (-) a+a 



(A2) 



(-) ata (a + a t )(-) QtQ = -(a + at) (A3) 

we obtain for Hh = \JHhU~ 1 

H H {t) \ 



H H = 



H H (-t) 



where 



H H (±t) = u> a j a - g(a + a 1 ) T t{-) ah 



(A4) 



(A5) 



In each block we have a purely phononic hamiltonian 
Hji(±t). The eigenvalues and eigenvectors can be de- 
termined by continued fraction solution for the resolvent 
between \m) and |n)-phonon states 



Gw„(w) = (m| 



1 



L0-H H (±t) 



Using 



1 



to — H u> — Hn lo — Hn u> — H 



(A6) 



(A7) 



with Hq = ujQa^a =F t(— ) a a and Hj = —g(a + a') we get 
the recursion 

-g^Gi(u>-moo)X ntP G^ n (uj) (A8) 
p 

where X U:P = (n\a + <v\p). This tri-diagonal recursion 
can be solved for the diagonal elements through a con- 
tinued fraction solution |E9(| 



Gt n (oj) 



1 



LO — ULOq ± t — S em — T, a bs 



(A9) 



where 

Z-'abs = 



ng 



LO + uo a =F t — 



(n - l)g 2 



lo + 2wo ± t 



(n-2)g 2 



LO + TILO0 + {=f) n t 

(A10) 



and 








(n + l).g 2 


LO - 


-LOQ^ft- 




{n + 2)g 2 




LO - 


... ±t (" + % 2 




lo — tiloq =p t — • • • 



and the property 



(All) 

At zero temperature the Green function of the two site 
Holstein model defined as 



10 



Gij(u>) = -t<0|Tc(t)c}(0)|0) 

can be expressed as 



Gi, 2 (w) = \(G+ (u) - Go i0 (w)) 



(A12) 



(A13) 



therefore Gq are the fc-space propagators whose poles 
determines the bonding and anti-bonding eigenvalues of 
Hh . The residues of the lowest energy pole of G^ n deter- 
mine the square of the projection of the phonon ground 
state on the |m)-state b^. Let us write Eq. (All) for 
n = in a recursive fashion 






pg 



uj - puj + {T) p t - 



(A14) 



j p+i 



where E^ is the continued fraction of Eq. ( All ) therefore 
the equation which gives the eigenvalue of the two site 
Holstein model is 

u±t-Bf = 0. (A15) 

By linearizing this recursion around the m-th solution 



E m of Eq. ( |A15| ) letting zf = d^/dE m we get 

-,2 



= «+l - 1) 



pg 



E m - ploo + {T) p t - E 



P+i 



(A16) 



. Finally the coefficient b^ is given by 



(A17) 



• \B) = \n)\i.),m = \n)\-i) 

. |C) = |M)lt-),|C") = |M)|.t) 



• \D) 



n)\i-),\D') = \n)\-i), 



where | ff) (| 4)) represent an up (down) spin state 
for the core spins and | |) (| J,)) are the same for the 
e g electrons. The Hamiltonian is invariant for flipping 
of all the spins so these are all the states we need. The 
states | A) and \D) have a FE character while the states 
\B) and \C) have AF character. The S = 3/2 subspace, 
spanned by the combinations of \A) and \A'), decouples 
from the other states even in the presence of e-ph phonon 
interaction 

If we consider the symmetric and antisymmetric com- 
binations 

|A±) = -L(|A) ± \A')) 
\B±) = -±=(\B)±\B>)) 

|G ± ) = -J=(|G)±|G')) 
\D ± ) = ±={\D)±\D')), 



(Bl) 



in the absence of e-ph interaction the symmetric and 
antisymmetric sectors are decoupled so that the part of 
the Hamiltonian matrix which pertains to Hund and an- 
tiferromagnetic interactions consists of 3 blocks: 



H- 



3/2 



Ji-Jh 



-t 




Ji-Jh 



+ t 



(B2) 



The (quantum) probability distribution function for 
the displacement operator can be determined using the 
harmonic oscillator wave functions \l/ n (X) as 



P(X) 



E 



[(b 



)*b + 



{b^)*b-} ** m (X)* n (X) (A18) 



APPENDIX B: THE EXACT SOLUTION OF THE 

TWO-SITE HOLSTEIN DOUBLE-EXCHANGE 

MODEL FOR S = 1/2 

Let us start form the case g = 0. We choose the fol- 
lowing basis set labelling the states according to the total 
spin Stot = S + S where S is the spin of the Mg 3+ ion 
and s that of the e g electron: We have two states in the 
S = 3/2 sector: 

• \A) = \n)\v),\A') = \n)\.i) 

and six states in the S = 1/2 sector: 



#1/2 



Ji + Jh 
2 L 





t 



2 

Ji — Jh 
4 

-is. 

2 





_ Jh 
2 

Jl+Jff 
4 



(B3) 



t 



the last block B\ii- can be obtained from Eq. 
(B3) by the substitution t — > — t. The e-ph ma- 
trix elements couple the subspaces (A + ,B + ,C + ,D + ) and 
(A~ ,B~ ,C~ ,D~). The subspace spanned by \A ) can be 
diagonalized independently having the same eigenvalues 
and eigenvectors of a 2 site Holstein model (see Appendix 
|A|). The hamiltonian matrix in the S\/-2 can be written 



_ / H 1/2+ +w la t a -gl(a t + a) 
^l/a- I _„iM + a ) H 1/2t . 4-^i«t, 



-gl(a) 



wola^a 



(B4) 



Here 1 is the 3x 3 unit matrix. We can diagonalize 
H1/2 in the phonon space by means of the unitary trans- 
formation 



i (i + (-r ta )i (-i + (-r ta )i 

2 I (-l + (-) ata )l (l + (-) ata )l 



(B5) 
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Using the property given in Eq. (A3) we get for -H1/2 = 

- / (uWa ~ 5(« f + o))l + Hf /a - t(-) a ' a A 

\ (woa'a- 

(B6) 

where we have splitted Hi/2,± — H"j 2 ±tA in the atomic 
(t = 0) part and in the hopping dependent term 



(B7) 




and -ffw n can be obtained from Eq. (B4) with t = 0. 



'1/2 



The hamiltonian Eq. (B6) can be diagonalized in each 
spin sector by an inddependent transformation which di- 
agonalizes -Hi/2,± for a given phonon number n. The 
diagonalization gives six eigenvalues Ei/2 t ±(n) f° r each 
phonon number. For n = and Jh S> t, J\ the lowest of 
them are 



(B8) 
(B9) 



By replacing t — ► £(— ) a a we are left with a purely 
phononic hamiltonian in the FE/AF sector 



FE 
1/2,- = 


4 


-< + f 


AF 
1/2,- ~ 


4 


£ Ji 
" 2 ~ Y 



H 1 / 2 — Uoa'a — g(a + a) 
H 1 , 2 = UJ o a 0, — g(a* + a) 



J_h_ 

4 

Jh_ 

4 



4 



2 l J 



■h 

2 



(BIO) 
(BH) 



The comparison between the energies of the FE and the 
AF phases leads to the condition (j|) for the magnetic 
transition. 
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